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A movement is underway to test the uniqueness of Earth. Sponsored primarily by NASA, 
it is enhsting talented researchers from many disciplines. It is conceiving new telescopes to 
discover and characterize other worlds like Earth around nearby stars and to obtain their 
spectra. The goal is to search for signs of biological activity and perhaps find other cradles 
of life. 

Most effort thus far has focused on the optics to make such observations feasible. Rel- 
atively little attention has been paid to science operations — the link between instrument 
and science. Because of the special challenges presented by extrasolar planets, science- 
operational issues may be limiting factors for the "Are We Alone?" (AWA) movement. 
Science-operational metrics can help compare the merits of direct and astrometric planet 
searches, and estimate the concatenated completeness of searching followed by spectroscopy. 
This completeness is the prime science metric of the AWA program. Therefore, the goals 
of this white paper are to present representative calculations involving science-operational 
metrics, and to promote a science-operational perspective. We urge the Survey Committee 
to allow this perspective and such metrics to inform its plan for the future of AWA. 

Extreme targets 

Because of three features, an Earth twin (ET) — the AWA benchmark adopted for the 
calculations presented here — is the most extreme target in astronomy: low information rate, 
high variability in position and brightness, and the huge, proximate, unwanted signal of 
the star. To illustrate, consider a space telescope with aperture Z) = 16 m, the largest 
currently under serious discussion. With efficiency e = 0.2, and resolving power i? = 10 in 
the / band, this telescope would detect 0.1 photons per second from an ET at quadrature 
at 20-pc distance. The apparent separation of the stellar and planetary images would be 
0.05 arcsec or 5\/D. Assuming an inner working angle IWA = 3X/D — meaning the radius 
of the photometrically inaccessible region immediately around the star — and assuming a 
maximum feasible difference in magnitudes between the planet and star Amago = 26, set by 
wavefront instability for an internal coronagraph, or by the starshade accuracy and alignment 
in the case of an external occultor, such an ET on a random orbit would be detectable 66% 
of r^mdom search epochsQ The median time tor a detected ET on a partially obscured 



^See calculations of completeness and discussion of AmagQ in Brown, R. A. 2005, "Single- Visit Photo- 
metric and Obscurational Completeness," ApJ, 624, 1010-1024. 
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orbit to become undetectable again — too faint or too close to the star — is two months after 
detection. Assuming the optics suppress the diffraction wings by factor 10"^'' from the 
central intensity of the stellar image, and assuming sharpness \E' = 0.02, some 10 photons of 
scattered starlight and 1 photon of zodiacal light — assuming 1.0 local zodi and 2.0 exozodi — 
would be counted with each planetary photon during the exposure time T = 146,000 sec that 
would be required to obtain a spectrum with photometric signal-to-noise ratio SNRpi^o = 10 
at the oxygen A band (760 nm), as required to measure 20% of Earth's column density of 
O2 with 99% confidence^] 

Turning to the astrometric signal, at 20 pc the amplitude of the astrometric wobble of 
the Sun due to the orbiting ET would be a = 0.15 /las or 0.00064 solar radius, a varies 
directly with the square root of the stellar luminosity and inversely with the stellar mass. 
An astrometric data set comprises measurements of the form (tj, Tj, Xj, yi, ai): the x- 
y position at epochs ti spread over more than one planetary period, where the exposure 
time is 0.5 for each of the two quasi-orthogonal directions, and where cxj = cro(ro/rj)^/^ 
is the positional uncertainty in either direction. The astrometric signal-to-noise ratio is 
SNRasi = a/a, where 

1/2 /v^^ Y^^^ 1/2 -1/2 



where r is the total exposure time. This noise reduction with r and A^ is valid only down 
to the noise floor, or for a > cifloor, where cxfloor is determined by systematic effects. For the 
Space Interferometry Mission {SIM Lite), a space astrometer currently under study, we use 
the instrumental parameters are do = 1.41 /xas, tq = 2200 sec, and cXfioor = 0.035 yuas|f|000 
Assuming periodogram (power spectrum) analysis of the data set to find weak, periodic 
signals, the search completeness C{SNRast, fap) depends on SNRi^t and the accepted false- 
alarm probability fap; for example, C = 0.5 for SNR^t ~ 6. The sensitivity limit for 
detection probability C = 0.5 is therefore amin = 6 x cTfjoor = 0.21 /ias, which would be 
achieved in r = Tmax = 3.6 x 10^ sec. For example, an ET at 14.3 pc would be just detectable 
by SIM Lite under these assumptions, because the required a = 0.21/6 = 0.035 fias, which 
is equal to the noise floor. 



^For details of the photometric calculation, see HandCheckPhotonietry_SunEarth.doc at 
http: / /homepage. mac.com/planetfinder /| 

■^Unwin, S. C, et al. 2008, "Taking the Measure of the Universe: Precision Astrometry with SIM 
PlanetQuest," PASP, 120, 38. 

''Catanzarite, J., Shao, M., Tanner, A., Unwin, S., & Yu, J. 2006, "Astrometric detection of terrestrial 
planets in the habitable zones of nearby stars with SIM PlanetQuest," PASP, 118, 1319-1339. 

^Traub, W. A., et al. 2009, in Extrasolar Planets in Multi-Body Systems: Theory and Observations, Ed. 
Krzysztof Gozdziewski, European Astronomical Society Publication Series, in preparation. 

^Brown, R. A. 2009, "On the Completeness of Reflex Astrometry on Extrasolar Planets near the Sensi- 
tivity Limit," submitted to ApJ; |http: //arxiv.org/pdf/0901 .489"77 
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Fig. 1. — AWA stars. Maximum apparent separations and astro metric amplitudes for an 
Earth twin (same size and temperature as Earth) around the 7322 stars in a SIM Lite 
project star hst. Colored stars (505): resolved at 3 A/D (blue line) and offering exposure time 
T < 500,000 sec for achieving SNRp = 10 at i? = 70 for the A band of O2 at A = 760 nm using 
a 16-m telescope. Red stars (117): ET astrometric wobble a > a^in = 0.21 /xas (red line), 
the sensitivity limit of SIM Lite (search completeness C = 0.5 and false-alarm probability 
fap = 0.01). Circled red stars (33): targets in the SIM Lite program that maximize ET 
discoveries (summed C) for 1 year of exposure time (100% duty cycle). Yellow: Sun-Earth 
at 20 pc. 

The mother of all problems 

The mother of all science-operational problems for AWA is the density of stars in the 
solar neighborhood. If it were higher, Hubble {D = 2.4 m) might have obtained the first 
spectrum of an ET — if it were capable of 10"^*^ contrast, of course. 

Figure 1 shows the 7332 stars in a SIM Lite project star list for planet searchinglll For 



^Thanks to Joe Catanzarite for developing and sharing this star list, which includes all the dwarfs and 
subgiants within 100 pc in the NStED data base (http://nsted.ipac.caltech.edu/) with stars cut by these 
criteria: (1) known from the Washington Double Star Catalog to have a companion closer than 35 mas, 
(2) known to have a companion between 35 mas and 1.5 arcsec with V magnitude difference <1, (3) lumi- 
nosity > 25, (4) fainter than V = 9, (5) chromospheric activity index Rhk > 4.3, (6) known to have age 
<0.3 Gyr, or (7) with orbital period at Earth's equilibrium temperature >4 years. 
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this list, a telescope with D = IQ m resolves 505 ETs at some time in their orbits and 
achieves the required SNRp = 10 for the O2 spectrum in less than T^ax = 500,000 sec, an 
arbitrary, "stretch" exposure time. ("ET" now means the same radius and albedo as Earth, 
but with semimajor axis equal to where L is the stellar luminosity, which ensures the 
same equilibrium temperature.) 

The exposure time for the O2 spectrum for the hypothetical Earth-Sun system at 20 pc 
is close to the median value for the 505 stars: Tmedian = 150,000 sec. 

Some 117 of the 505 stars suitable for O2 spectroscopy with a 16-m telescope have ETs 
detectable by SIM Lite according to the sole criterion a > amin = 0.21 //as. 

Finding their addresses 

For ETs, the two main searching options under consideration today are direct imaging 
and astrometry. Radial velocity works for larger planetary masses, but probably not ETs. 
For example, the 505 stars just discussed have a median mass of 0.97 and range from 0.4 to 
2.4 solar masses. The radial-velocity amplitude for the Earth-Sun is 0.089 m/sec, assuming 
an edge-on orbit. For other stars, the ET radial- velocity amplitude varies inversely with both 
the square root of the mass and the fourth root of the luminosity. For the 505 stars, the 
median radial-velocity amplitude is 0.094 and ranges from 0.03 to 0.3 m/sec — substantially 
beyond current radial-velocity performance limits. Occultation photometry won't work be- 
cause the alignments are highly improbable for the number of nearby stars. 

Discovery by direct detection was studied intensively for TPF-Cj^ It tends to be ineffi- 
cient — low completeness per search observation due to the highly obscured orbits, and long 
exposure times — but still feasible. Certainly, any telescope that could obtain the O2 A-band 
spectrum of an ET at i? = 70 could detect the planet by direct imaging in a shorter time 
at i? = 3-10. Nevertheless, a search calls for starlight suppression over a large field of view, 
whereas spectroscopy requires it only at the planetary position — if the position is known. 
Furthermore, direct detection has other problems that have not been solved. For example, 
it cannot easily determine the identity of a faint image, which could be an ET, a speckle, a 
background source, a cluster of exo-zodiacal light, or a companion of only indirect interest for 
AWA, such as a Neptune or Jupiter. The discovery image — and indeed all the data obtained 
in the first observing season — would provide neither the mass nor the orbit. (Solar avoidance 
angles are typically >90°, which means observing seasons <6 months; orbital coverage is 
inadequate in the first season.) First-season images could disambiguate the only a fraction 
of sources, depending on a star's particular combination of parallactic ellipse and proper 
motion. Another debilitating consequence of not determining the orbit is that recovery of 



^Brown, R. A., Hunyadi, S. L., & Shaklan, S. B. 2006, "A DRM for TPF-C: A Design Reference Mission 
for the Coronagraphic Terrestrial Planet Finder;" 
http://planetquest.jpl.nasa.gov/TPF/tpfcDocs/Brown.pdf/ 
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the planet in future observing seasons — predicting its position for planning and scheduling 
purposes — becomes a low-probability, hit-or-miss proposition!^ 

To further investigate astrometric discovery, we conducted a reverse auction of exposure 
time r to maximize the number of planets discovered (total completeness) as a function of 
exposure time on SIM Lite (assuming 100% duty cycle). First, = Xmax was assigned 
to each of the 505 stars. Next, SNR^st and the completenesses Ck were computed using 
the "discovery" curve in Fig. 1 of Ref. 6. Next, those steps were repeated for minus 
2200 sec (a quantum of exposure time), and the star was identified for which the reduction in 
completeness would be smallest. For that star, the reduction in and associated reduction in 
Ck were implemented accordingly. This process was repeated many times, always removing 
the smallest amount of completeness per 2200 sec of exposure time, until the point was 
reached where Sta; = 0. The auction record is shown in Fig. 2. At Sr^ = 1 year, 33 stars 
were still being observed. For these stars, the minimum, median, and maximum values of 
SNRa were (5.8, 7.2, 11.1), with corresponding values of values of C: (0.44, 0.74, 0.99). 

The discovery of a planet — yes or no — is a Bernoulli random variable with probabil- 
ity C. The expectation value of the sum of Bernoulli random variables is the sum of their 
probabilities. By summing the values of Ck, we find the expectation value of the number of 
ETs discovered. For Sxfc = 1 year, SCfc = 24.6. In this case, if the occurrence rate of ETs is 
rj = 0.05, the most probable result would be one ET discovered but a 30% chance of zero. 

For Sr^ = 12 years, all a = (Tuiin, all 117 stars are searched to their limits, and SC^ = 

94.7. 

Astrometry's unique forte is determining the masses and orbits of the ETs discovered. 
However, the least-squares estimators of Keplerian parameters, particularly for mass and 
eccentricity, are biased in the low SNRg^t regime, which is the domain of ETs for SIM LiteS 
Preliminary indications are that, as a result of this bias, the completeness for mass estimates 
as good as ±25% would be below 0.7, and below 0.3 for ±10% accuracy, for typical ETs 
observed by SIM Lite. The biased estimation of orbital elements also reduces the accuracy 
of positional estimates.^ 

Discussions 

The above calculations are science operational. They estimate the science produced 
by operating the hardware. Because science is the ultimate purpose of the hardware, such 
calculations are essential for defining, understanding, and prioritizing exoplanetary projects. 
This is particularly important for the AWA movement, both because of the dynamic, low- 
SNR targets, and because complementary, concatenated missions may be involved. 



^Brown, R. A., Shaklan, S. B., & Hunyadi, S. L. 2006, "Science Operations and the Minimum 
Scale of TPF," JPL 07-02 7/07, Proceedings of TPF-C Workshop 9/28-29/2007. (SciOpsTPF.pdf at 
|http : / /homepage . mac . com/planetfinder / , ) 

^°For astrometry, see Ref. 6; for radial velocity, see: Shen, Y. & Turner, E. L. 2008, "On the Eccentricity 
Distribution of Exoplanets from Radial Velocity Surveys," ApJ, 685, 553. 
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Fig. 2. — Auction results. Red: SIM Lite total search completeness SC^ and stars in play 
(still in the auction), versus total exposure time in an observing program that maximizes 
SCfe. For 1 year total exposure time (100% duty cycle, 33 stars are in play, and summed 
SCfc = 24.6, which is the expectation value of the number of ETs discovered, assuming all 
stars have one. In 12 years of exposure time, all 117 stars with ETs detectable by SIM Lite 
are in play, and the maximum summed completeness is achieved, SCfc = 94.7. 

Because of unresolved science-operational issues, it is challenging to outline a coherent, 
low-risk AWA mission series, progressing systematically from discovery of ETs to atmospheric 
spectroscopy in search of biomarkers. The results of Kepler will certainly help by clarifying 
what total search completeness is required to ensure an adequate sample of ETs. Even when 
the occurrence rate of ETs is better known, however, we still need deliberate processes for 
discovering, identifying, and obtaining the ephemerides of nearby ETs. 

That is the top-down approach. The default, bottom- up approach would have the AWA 
movement not defining single-purpose missions, but proceeding opportunistically, competing 
for instrumental and observational opportunities on missions that are justified primarily by 
other science priorities. For example, the TPF-C studies illustrate that high-performance 
coronagraphy is compatible with wide-field imaging — not only for optics and instruments on 
an off-axis, unobscured telescope, but also in terms of science operations!^ 



^^R. A. Brown et al. 2006, "Final Report of an Instrument Concept Study for a Wide-Field Camera for 
TPF-C' See WFC ICS Final Report.pdf at http: / /homepage.mac.com/planetfinder/ . 
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Someday, Hubble^s successor will be built to rekindle research along the long frontiers of 
science for which Hubble defines the present state of knowledge. A coronagraph will surely 
be on board — this time with realistic hopes for the goal of "direct imaging . . . for planetary 
companions of nearby stars. "0 To illustratate the relative power of direct imaging, we note 
that all 33 stars circled in Fig. 1 are bright and close. Their ETs would be resolved with 
an IWA = 3X/D for an 8-m telescope, which would achieve SNRpi^o = 10 at i? = 5 in the 
I band with a median exposure time of only T = 3.4 x 10^ sec. For 27 of the 33 stars, a 4-m 
telescope would do the same with a median exposure time T = 5.6 x 10^ sec. 

Similarly, in the case of space astrometry, astrophysical science other than exoplanetary 
was strongly endorsed by two previous decadal surveys, and it remains the SIM Lite mission's 
primary justification. If SIM Lite finds ETs, it will prove the existence theorem and produce 
at least rough estimates of mass and orbit. This will be of great interest. As illustrated in 
Fig. 2, more observing time would reach more stars: 12 years of SIM Lite exposure time, or 
12 SIM Lites each exposing for 1 year, would discover all the ETs available to SIM Lite in 
this star sample. 

Final remarks 

The fact we can even contemplate such challenging observations attests to NASA's 
progressive mastery of ever-broader wavefronts of astronomical light. This great achievement 
is evidenced by the success of the Hubble Space Telescope, as well as by plans for the James 
Webb Space Telescope, as well as designs for the SIM Lite and the coronagraphic Terrestrial 
Planet Finder (TPF-C). This mastery, on which any optimistic view of the long-term future 
of astronomy must largely be based, will be counted among the finest achievements of the 
space program. 

In his memoir, Riccardo Giacconi offers a comment pertinent to AWA, on "the greatly 
increased program of research on extrasolar planets, stemming mainly from the personal 
interest of Daniel S. Goldin (a manager and engineer), NASA's administrator from 1992 to 
2005. This increased effort was welcomed by a minority of the astronomical community, but 
happened without prioritization with respect to other fields in astronomy of much greater 
astrophysical interest."!^ Hopefully, in the current decadal study, the AWA movement will be 
reconciled with other priorities, based on the AWA science that could actually be produced 
by operating the hardware, as estimated by science-operational metrics. 

I thank Neill Reid, Marc Postman, Roeland van der Marel, and Remi Soummer for their 
comments, and Sharon Toolan for helping to prepare this paper. It was originally written 
as a science white paper for Astro2010. 



^^A listed science objective in NASA's "Announcement of Opportunity for Space Telescope," March 1977. 
^•^Giacconi, R. 2008, Secrets of the Hoary Deep: A Personal History of Modern Astronomy, (Baltimore: 
The Johns Hopkins University Press). 



